Creep of a fracture line in paper peeling 
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The slow motion of a crack line is studied via an experiment in which sheets of paper are split into 
two halves in a "peel-in-nip" (PIN) geometry under a constant load, in creep. The velocity-force 
relation is exponential. The dynamics of the fracture line exhibits intermittency, or avalanches, 
which are studied using acoustic emission. The energy statistics is a power-law, with the exponent 
(3 ~ 1.8 ± 0.1. Both the waiting times between subsequent events and the displacement of the 
fracture line imply complicated stick-slip dynamics. We discuss the correspondence to tensile PIN 
tests and other similar experiments on in-plane fracture and the theory of creep for elastic manifolds. 
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The deformation and fracture of materials is a fasci- 
nating topic as one can explore the physics even without 
sophisticated experiments [l|. A piece of paper suffices 
to give ample evidence for the presence of phenomena 
that need a statistical description. One can tear a sheet 
or crumble it, to observe that the response is "intermit- 
tent" and not simply "smooth". [1 H 0, [1] . 

The physics of fracture in a material as paper is gov- 
erned by two basic ingredients. The structure and the mi- 
croscopic material properties are inhomogeneous, while 
the stress fields follow the laws of elasticity. Typical sta- 
tistical signatures that show features that emerge from 
their interaction are the acoustic emission during a de- 
formation test and the post-failure properties of fracture 
surfaces. These are commonly found to be described by 
power-laws, as regards the energies of acoustic events 
("earthquakes"), the intervals between subsequent events 
( "Omori's law" ) [3, 0] , and the same is true for the geome- 
try of cracks where ample evidence points towards a self- 
affine fractal scaling of the surface fluctuations around 
the mean, which is found in a variety of cases (loading 
conditions, materials and so forth) 0, HQ- 

The most simple case of fracture - which we also will 
study here - is the passage of a crack line through a sam- 
ple, when its movement is constrained on a plane. In this 
case, the mathematical description is given by a crack po- 
sition h(x,t), where x is along the average projection of 
the crack. The average motion is described by h = vt 
with v the crack velocity. In the situation at hand one 
has three ingredients: a disordered environment which 
poses obstacles to the line motion, a driving force K ext 
(stress intensity factor in fracture mechanics language), 
and the self-coupling of the fluctuations in h through a 
long-range elastic kernel (lp| . which is expected to scale 
as 1/x. 

To describe the line's physics one needs the language 
of statistical mechanics. One finds a phase diagram for 
v{K ext ): an immobile crack begins to move at a critical 
value K c of K e ff such that for K > K c v > 0. This 
transition between a "mobile line" and a pinned line is 
commonly called in non-equilibrium statistical mechan- 



ics a "depinning transition". Close to the critical point 
K c the line geometry is a self-afhne fractal with a rough- 
ness exponent £. This is an example of a wide class of 
similar problems, which range from fire fronts in com- 
bustion to flow fronts in porous media to domain walls 
in magnets. The planar crack ll| or the contact line 
on a substrate [13] problem has been studied both theo- 
retically via renormalization group calculations and nu- 
merical simulations, and via experiments. There is a dis- 
crepancy between these, in that the roughness exponents 
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are (theory ~ 0.39 and C, exp t ~ 0.6 [lc 
dimensional experiments the in-plane roughness has re- 
cently shown signatures of ( expt ~ (theory UM- Imaging 
experiments have shown the existence of avalanches, of 
localized intermittent advances of the crack front with an 
avalanche size s distribution P(s) ~ s -i-6--i.7 [TJ| 

Our work considers the dynamics of such a crack front 
during creep, which is done by peeling paper sheets in 
a geometry illustrated in Fig. [T] (see also [Hj]). The 
creep of elastic lines (or manifolds or domain walls) 
is important since when K e fj_ < K r thermal fluctua- 
tions take over for T > [HUl O- The fluctua- 
tions nucleate "avalanches" similarly to what happens 
in zero-temperature depinning in the vicinity of K c . The 
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avalanches induce a finite velocity v, 
ment of a crack front might be fundamentally different 
from say a domain wall or a contact line, which can fluc- 
tuate back and forth around a metastable state: it is 
easy to see that the crack may in many cases only grow 
forward. 

The form of v creep (K e ff,T) depends on the "energy 
landscape" since the current understanding is that c reep 
takes place via nucleation events over energy barriers [l7| . 
These barriers are related to the roughness exponent ( 
and to its origins. The important physics is summarized 
with the creep formula 



exp-(C/m AI ) 



(1) 



which states that the creep velocity is related to the driv- 
ing force, m (as for mass, see Fig. [T] again),. There is a 
creep exponent /i, which depends on the interactions and 
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dimension of the moving object (a line). One expects the 
scaling 



fi = e/v 



2C 



(2) 



in d dimensions, (d — 1 for a line, d = 2 for a domain 
wall in say a bulk magnet). 9, v, and C denote the energy 
fluctuation, correlation length, and roughness exponents 
relevant for the problem, and exponent relations are usu- 
ally postulated between these three. They exponents all 
depend on a, the fc-space decay exponent of the elas- 
tic kernel. For long range elasticity, one would assume 
a = 1 whereas for the so-called local case a = 2 is ex- 
pected. Expressed as above, the question of the value 
of fi boils down to the values of £ and a - what is the 
effective roughness at hand, and what is the elasticity of 
the line like? 

The fundamental formula of Eq. @ has been con- 
firmed in the particular case of 1+1-dimcnsional do- 
main walls [20( | . and further experimental studies exist 
PH . Our main results on fracture line creep concern the 
velocity- force- relation v(m), and on the picture of the dy- 
namics that ensues, in particular from Acoustic Emission 
(AE) data of the avalanches and their dynamics. We find 
an exponential v(m) and discuss its interpretations. The 
dynamics shows signatures of intriguing but weak corre- 
lations, and we discuss the similarities and differences to 
the tensile case. 

In Figure [1] we show the apparatus • The failure 
line can be located along ridge, in center of the the Y- 
shaped construction formed by the unpeeled part of the 
sheet (below) and the two parts separated by the ad- 
vancing line. Diagnostics consist of an Omron Z4D-F04 
laser distance sensor for the displacement, and a stan- 
dard plate- like piezoelectric sensor [22| ■ It is attached to 
the setup inside one of the rolls visible in Fig. [1] and the 
signal is filtered and amplified using standard techniques. 
The data acquisition card has four channels at 312.5 kHz 
per channel. We finally threshold the AE data. The dis- 
placement data is as expected highly correlated with the 
corresponding AE, but the latter turns out to be much 
less noise- free and thus convenient to study. For paper, 
we use perfectly standard copy paper, with an areal mass 
or basis weight of 80 g/m 2 . Industrial paper has two 
principal directions, called the "Cross" and "Machine" 
Directions (MD/CD) (23[. The deformation characteris- 
tics are much more ductile in CD than in MD, but the 
fracture stress is higher in MD. We tested a number of 
samples for both directions, with strips of width 30 mm. 
The weight used for the creep ranges from 380 g (CD) to 
533 g (MD), and the resulting data has upto to tens of 
thousands of AE events (avalanches) per test (Fig. [3]). It 
is unforunately not possible to infer the critical depinning 
m c > m use d- The mechanical (and creep) properties of 
paper depend on the temperature and humidity. In our 
setup both remain at constant levels during an experi- 



ment, and the typical pair of values for the environment 
is 49 rH and 27 °C (though we also have data for other 
combinations). 

Our main result is given in Figure [H where we show 
the v(m) vs. 1/m characteristics. There are four main 
datasets depicted. These differ slightly in the typical 
temperature and humidity (for set 2 44 % rH instead of 
about 49 % rH for the others). These all imply an ex- 
ponential behavior, and by fits to the data sets we can 
infer a creep exponent v = 1.0 ± 0.1. This lends itself to 
two different interpretations: either the creep is in the Id 
random field (RF) domain wall universality class, since 
with a = 2 v — 1 implies C — 1, the roughness expo- 
nent of the RF universality class at very small external 
fields/forces. This assumes local elasticity being domi- 
nant. If we would then assume that £ takes the value 
C ~ 0.6 . . . 0.7, of the avalanche regime, this would imply 
screened interactions with a ~ 1.4. The thirdr possibil- 
ity is to use non-local line elasticity, with a = 1. Then 
Eq. [2] indicates that C ~ 1/3- This is exactly the equi- 
librium exponent of lines with long-range elasticity [261 ] . 
The velocity in the case of "paper" is influenced by the 
humidity: this is clear in our case. Meanwhile, note that 
the temperature variation is insignificant here. 

The avalanche behavior is illustrated by Fig. [3^, by 
the way of the AE event energy distribution P(E). We 
present two kinds of avalanche data: integrated over 1 ms 
windows (Ei where i is an integer time index for CD) and 
extracted from single avalanches. The former obviously 
sums over all the avalanches during the 1 ms duration 
(if more than one are present). The data agrees rather 
well regardless of the mass m and the ductility (CD, MD) 
with the scaling P(E) ~ E^ 13 with (3 = 1.7 ± 0.1. This 
is close to the one observed in the corresponding tensile 
experiment (depicted in the Fig. 0k) [Hj]. The simplest 
interpretation would be that once an avalanche is created 
through a thermal fluctuation, it follows a deterministic 
course. This is similar to the predictions of recent the- 
ories, which indicate the presence of a small nucleation 
scale, and that the avalanches should be as at the depin- 
ning critical point [H, [H, E3] . 

We have studied the temporal statistics and correla- 
tions using the AE timeseries, in particular the windowed 
one, and the direct displacement vs. time -signal (note 
that the AE has better accuracy in the time domain). 
The waiting times r, the silent intervals between either 
two avalanches or two windows with Ei > are shown 
in Fig. UJd. It is interesting to note how the distributions 
P(t) ~ r a change with the applied force. For large driv- 
ing forces it appears that there are (perhaps) two power- 
law like regimes: one with an exponent a — 1.3, and more 
clearly one for large r with a ~ 2 which is also found for 
m small. The general form for m large resembles also 
that of similar tensile tests where a ~ 1 is found. In the 
tensile case, there is a typical time/lengthscale (arising 
from paper structure) visible in tensile crackline tests, 
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which might here also be related to the change in the 
slope of P(t). 

The dynamics of the line exhibits stick-slip: the in- 
tegrated velocity depends on the window length under 
consideration. This is already evident from the P(E) 
and extends to longer timescales than what P(E) im- 
plies via durations of events. A "stick-slip" exponent can 
be extracted also from the displacement data as well as 
integrated from the AE data. It appears that a relation 
P(Ah) ~ A/i~ 17 arises. This implies the presence of 
correlations, which we next study by the energy release 
rate T: the duration of the active time T as measured 
as subsequent windows where Ei > 0. It can be seen as 
illustrated in Figure [5] that P(T) ~ T~ 2 - 7 , so that on a 
millisecond scale (much slower than avalanche durations, 
but much faster than the inverse line velocity) there are 
clear correlations. These are however weak in the respect 
that one can not see in here or the AE data signatures of 
"aftershocks" or "precursors" familiar from earthquakes 
or from the AE activity of dislocation systems [25]. One 
may also see a correlation between the energy released in 
the active period vs. its duration, such that E to t ~ T 25 
at least for m small. Thus it appears that the tempo- 
ral dynamics can be described by a temporal self-affinc 
process. 

To conclude, we have studied the creep dynamics 
of an elastic line, or a fracture front in peeling paper 
sheets. The particular features of our case are the disor- 
der present in usual paper and the variation in material 
properties (ductile/brittle). The main findings are four- 
fold. First, we have obtained an exponential velocity- 
force relation, which has three interpretations; we prefer 
the one which assumes non-local elasticity implying that 
an equilibrium £ ~ 1/3 governs the creep. Second, we 
observe intriguing similarities and differences in AE or 
avalanches to ordinary tensile (constant, small velocity) 
experiments. These, third, indicate the separation of de- 
terministic, zero-temperature dynamics from the nuclc- 
ation - as in creep indeed - of individual events. Fourth, 
the dynamics of the process exhibits correlations that 
would need a theoretical explanation. Our results clearly 
call for more theoretical effort to understand in-plane 
fracture fronts, and their creep properties. They also 
indicate the need for general studies of creep fracture as 
a statistical physics problem. Acknowledgments: The 
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FIG. 1: Schematic view of the experimental setup (rolls, pa- 
per, AE sensor, camera, weight), and two closer views of the 
geometry. 
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FIG. 2: An example of the acoustic or stick-slip activity (Ei > 
0) during one single creep test. The energy Ei is integrated 
after thresholding over 1 ms windows. 
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FIG. 3: The creep velocity vs. the inverse of the applied 
force or mass. The line indicates an exponential decay (v ~ 
expa/m). The four sets (circle,rectangle,diamond,triangle) 
differ such that the 2nd one was done under lower relative 
humidity, rH % 43). 




FIG. 4: a): The probability distributions of the acoustic 
events for various cases, b): pots of the event interval (wait- 
ing time) distributions, for different massses m. The data sets 
have been shifted for clarity. 
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FIG. 5: The P(T) of durations of active times T (for definition 
of T see text) The line has slope -2.7. 



